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Lattice fermions obeying the Ginsparg- Wilson relation do correctly represent the 
physical properties related to chirality. This can be achieved by local fermions, 
which involve an infinite number of couplings, however. For practical purposes, it 
is useful to first construct approximate Ginsparg- Wilson fermions within a short 
^sj ' range. We report on a successful construction in QCD at (3 = 6. The good quality 

of the approximation is observed from the spectrum, which is situated close to a 
Ginsparg- Wilson circle. These fermions also provide an excellent approximation 
to rotational symmetry and they are promising for a good scaling, since they arise 
from the perfect action framework. Their insertion into the overlap formula renders 
the Ginsparg- Wilson relation exact. It leads to an improved overlap fermion with 
a high level of locality. This insertion is statistically on safe grounds at /3 > 5.6. 



o 

£—n ' 1 Ginsparg- Wilson fermions 

In a slightly simplified form, the famous Nielsen-Ninomiya theorem El states 
that a local lattice fermion without species doublers cannot be chiral in the 
sense that the lattice Dirac operator D anti-commutes with 75 . Locality means 
here that the couplings in D decay at least exponentially in the separation 

CIh! between ip and ip. n Hence it is an obvious idea to break the chiral symmetry 

by an irrelevant term, so that it should be restored in the continuum limit. 
The simplest way to do so is to set ^j5{D,j^} equal to some local term of 
0(a 2 ), such as the Wilson term ^A, where A is a discretized Laplacian. [] 
However, this type of chiral symmetry breaking on the lattice is rather violent; 

jJJ ■ it causes quite some trouble such as additive mass renormalization, 0{a) scaling 

artifacts, renormalization of currents, mixing of matrix elements, etc. 

On the other hand, it turned out to be harmless to introduce a non- 
vanishing anti-commutator as 

\ l5 {D-\ ln } = R (1) 

where R is a local term with {i?,7s} ^ 0. The superiority of this relation can 
be understood intuitively from the fact that R doesn't shift the poles in D^ 1 . 

"Actually the proof in Ref. [nl] still holds for an even weaker form of locality. 

'Here a is the lattice spacing, but in general we will refer to a hypercubic lattice of unit 
spacing in Euclidean space. 
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In the form 

{D XtV , 75 } = 2{D l5 RD) x>y (2) 

it is known as the Ginsparg- Wilson relation (GWR) 0. 

1.1 Virtues 

Amazingly, it seems that all physical properties related to chirality are re- 
presented correctly by a lattice fermion obeying the GWR (a GW fermion). 
The mass and the vector current (as well as the flavor non-singlet axial vector 
current) are not renormalized, and weak matrix elements dp not mix El. More- 
over, the chiral anomalies Ho0 as well as global anomalies 0, and the soft pion 
theorems are reproduced correctly 0. Even the construction of chiral gauge 
theories on the lattice is feasible based on GW fermionsQ. 

For the understanding of these properties, it is a key observation that GW 
fermions have an exact — though lattice modified — chiral symmetry at finite 
lattice spacinga. 

It is instructive to consider the spectrum of a GW fermion. For simplicity, 
we assume D' — 75-D75 and R x ,y = 4~^x,y (m > 0), hence the GWR reads 
fj,(D + D+) = D^D. If we introduce the operator A = D - /1, the GWR 
simplifies further to 

A^A = n 2 . (3) 

Therefore we know that the spectrum of a GW Dirac operator is — with the 
above assumptions — always situated on a circle in the complex plane, with 
center and radius /x. This confirms the absence of additive mass renormaliza- 
tion, and it also rules out "exceptional configurations" a. Moreover, it provides 
a well-defined index (since the zero eigenvalues are exact), and together with 
the index theorem we obtain a sensible definition of the topological charge of 
a lattice configuration. 

1.2 Limitations 

After celebrating the impressive properties of GW fermions, we now have to 
address their limitations. One point is that the GWR guarantees a correct 
chiral behavior, but it does hardly imply anything about other properties, 
which are also essential for a formulation of lattice fermions, in particular the 
scaling behavior. 

A second point concerns locality: the relaxation of the condition {D, 75} = 
to the GWR allows the lattice fermion to be local in the sense that the cou- 
plings decay exponentially with the lattice distances. This is sufficient from 
a conceptual point of view, since there is a decay length of a finite number 



of lattice spacings, which ensures the right continuum limit. However, for ap- 
plications one would like to have even "ultralocality" , which means that the 
couplings drop to zero beyond a finite number of lattice spacings. Unfortu- 
nately, GW fermions cannot have this property, not even in the case of free 
fermions E3. For example, if we want to insert the free Wilson Dirac opera- 
tor and solve for R, then we obtain a pseudo-GW kernel R which decays as 
R x .y oc \x — y|~ 4 in d = 2, and like R x . y oc \x — y\~ 6 in d = 4. This is non-local, 
and therefore the Wilson fermion does not obey any GWR, not even in the free 
case (in the interacting case this is also clear from the mass renormalization). 
Of course, in practice one cannot work with couplings over infinite dis- 
tances. In a finite volume with certain boundary conditions, the GWR - 
with these boundary conditions implemented — can be solved, but this still 
requires the coupling of sites (and links) over all distances in this volume. 

1.3 Exact and approximate solutions 

Regarding the first limitation, there exists a class of lattice actions called "per- 
fect actions" which deserves its name by yielding a scaling identical to the 
continuum at any lattice spacing. At the same time, perfect actions solve the 
GWR El, but unfortunately their construction is about as difficult as solving di- 
rectly the model under consideration, since it requires a functional integration 
extrapolated to the continuum. 

The construction of "classically perfect actions" liil is much easier, though 
still difficult. They also solve the GWR LI, and their scaling is still excellent. 
However, a successful construction apd application for interacting fermions 
could only be achieved in d = 2 so farl!3. From the second limitation we know 
that such actions need some truncation. Here we truncate to a "hypercubc 
fermion", HF (with couplings not only to nearest neighbor sites, but to all 
es inside a unit hypercube), which is actually applicable in QCD simulations 
, We first consider the truncation for the free fermion. There the perfect 
lattice Dirac operator can be constructed, and the term R occurs in the block 
variable transformation. Locality is optimal — thatJs, the exponential decay 
of the couplings is optimally fast — for R x , y — \^x,y^ which we denote as the 
"standard GW kernel" . Hence the spectrum of that perfect fermion is situated 
on a unit circle with center 1. We now truncate by evaluating the perfect 
couplings in_a small volume of 3 4 sites, and then we use the same couplings in 
any volumelij. We obtain a lattice Dirac operator of the form 

D(x, y) = p^x - y)7 M + X(x - y) , (4) 

where the support of p^ix — y), X(x — y) is restricted to \x v — y v \ < 1, (y = 
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1...4). (These couplings are given in Ref. JTa], Table 1.) Of course, in a larger 
volume it is not exactly perfect any more, and the GWR is violated a little. To 
probe this truncation effect, we look at the spectrum of the hypercube fermion 
on a 20 4 lattice. Fig. [y shows that it is indeed very close to a GW circle. Also 



0.5 



-0.5 



„sr 



<*•*•"*•" 



*****«, 




0.5 



1.5 
Re 



2.5 



Figure 1: The spectrum of a truncated perfect, free HF on a 20 4 lattice (plotted in (D). 

the scaling artifacts in this truncated perfect free fermion are very small, as we 

Se_from the fermion dispersion relation and thermodynamic scaling quantities 

The same construction can be done in d = 2, .-but there we prefer to use 
a "scaling-optimal hypercube fermion" (SO-HF) 113, which is similar to the 
truncated perfect one, but still a bit improved (with respect to scaling) . Again 
the spectrum for the free fermion is close to a GW unit circle, and its scaling 
is excellent. To simulate this lattice fermion in the 2-flavor Schwinger model, 
we performed a "minimal gauging" by hand: we attached the free couplings 
to the shortest lattice paths only, in equal parts where several shortest paths 
exist. We also added a clover term with coefficient 1, and used the standard 
plaquette gauge action. Of course, this simple gauging brings in additional 
artifacts and a further deviation from the GWR. Hence the eigenvalues spread 
more and more around the unit circle as (3 decreases, but for instance at (3 — 2 
the circle is still approximated more or less, see Fig. 0. 

As a scaling test we measured the dispersion relations of the two meson- 
type states, a massless triplet and a massive singlet £3, which we denote as it 
and rj. Indeed, the scaling is dramatically improved over the Wilson fermion 
(at critical hopping parameter), see Fig. [|. Amazingly, this SO-HF reaches a 
similar scaling quality as the classically perfect action, which was truncated 
only very mildly to 123 independent couplings per site Ii3. In contrast, the 
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Figure 2: The spectra of the 2d scaling optimal hypercube fermion (SO-HF) for typical 
configurations at strong resp. weak coupling, approximating a GW unit circle. 

SO-HF involves only 6 independent couplings per site, hence this approach 
has the potential to be extended and applied in d = 4. 
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Figure 3: The "meson" 
Schwinger model at f$ = 6 
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Also the rotational symmetry is approximated very well for 



ir such HFs. 
i[e tested this for free and for interacting fermions in d = 2 IIZI and d = 4 
E3. As an example, we show the "speed of light" for the minimally gauged 
truncated perfect HF in QCD at (3 = 5, and compare it to the Wilson fermion 
in Fig. H. What is not visible from that figure, however, is the dramatic mass 
renormalization: what is supposed to be the pion mass amounts to M = 3.0. 
Also in the Schwinger model this problem is serious, as Figs. @ and H show: 
at = 6 we obtain a 7r mass of 0.13. This is the one unpleasant feature of 
the otherwise successful HFs with minimal gauging. In the next section we 
discuss a possibility to eliminate this effect completely. In Sec. 4 we are going 



to consider further methods to approach the chiral limit. 

2.0 | ■ , . , . , 



continuum 
• Hypercube 
i Wilson 



5 

z 



1.0 



1.0 

p 



Figure 4: The "speed of light" c = p/VE 2 - M 2 (M = "pion mass") in QCD at f3 = 5. E3 

2 Improved overlap fermions 

Let us start from some lattice Dirac operator Dq, which obeys the conditions 
of the Nielsen-Ninomiya theorem (local, no doublers etc.). We assume again 
Dq = 75-D0751 an d we recall that the GWR for R x y = -^-5 X:V is equivalent to 

A' Ao = /i 2 , for Aq = Do — (J,, (/x > 0). In general this will not hold, of course, 
but we can simply enforce it by the overlap formula 

' - / '" 1 " ' D ov = A ov +^. (5) 
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D ov represents a GW fermion. A prototype was proposed by H. Neubergerliia, 
starting from the Wilson fermion Dq = Dw and \x = 1, 

D Ne = 1 + ; W ■ A w = Dw - 1 . (6) 

A W A W 

In the free case and in a smooth gauge background, Djv e is free of doublers 
and local. The latter has heeiLestablished analytically assuming a constraint 
on each plaquette variable o'Eil, and numerically for QCD at /3 = 6 E3. For 
first simulations in quenched QCD, see Refs. [ p2|]23| ]. 

However, this is just one example iri a large class of overlap fermions, 
which is obtained by varying Dq and /i. L3 In particular, if it happens that Dq 
represents a GW fermions already, then it is just reproduced under the overlap 
formula, D ov — D (for a fixed GW kernel R). Therefore, also perfect and 



classically perfect fermions are special types of overlap fermions. In practice 
we do not have an exact GW operator Dq at hand, but we can construct an 
approximation, as discussed in Sees. 1.3 and 4. Then the square root will keep 
close to the constant /i, and D ov m Do- We may say that the overlap formula 
provides a "GW correction" of Do; m particular it removes the additive mass 
renormalization. 

Based on this property, we suggest the following concept a: start from a 
truncated perfect fermion (or something similar), which scales well and approx- 
imates the GWR, and insert it as Dq into the overlap formula. The resulting 
D ov has exact GW chirality (in particular we are in the chiral limit) and a 
good scaling behavior can also be expected, due to D ov rs Dq. The latter also 
suggests that an approximate rotational symmetry of Do is essentially inher- 
ited by D ov , and that D ov is very local in the sense of a fast exponential decay: 
remember that Dq is ultralocal, and — due to the modest modification — the 
long-range couplings will be turned on just a little bit in D ov . H 

To summarize, this concept aims at combining all desirable properties of 
a lattice fermion formulation. It has been tested comprehensively in the two- 
flavor Schwinger model, and the above expectations are observed to hold in a 
very satisfactory way. In rest of this Section we summarize our 2d results from 
Ref. [0]. 

We first checked the scaling quality of the free overlap fermions (fermionic 
dispersion, thermodynamic scaling ratios) and it turns out that the improve- 
ment of the SO-HF over the Wilson fermion persists under the overlap formula, 
i.e. the overlap SO-HF is strongly improved over the Neuberger fermion. By 
comparing the spectrum of a fixed configuration for the SO-HF and for the 
overlap SO-HF, we could then literally see that the alteration due to the over- 
lap formula is small; the eigenvalues are moved almost radially onto the unit 
circle. Next we tested the scaling in the interacting case, and the drastic 
improvement of the overlap SO-HF over the Neuberger fermion is confirmed 
again, see Fig. 0. We have used standard operators in both cases, hence the 
comparison is fair. Of course, both types of overlap fermions may improve if 
one consistently improves the operators, which can be_tedipus, however. Then 
one expects also the Mauberger fermion to scale better cSEij (in particular 0(a) 
artifacts are excludedEj) , but our results imply that improved overlap fermions 
can be used very successfully even with the simple standard operators (as it 
was also observed for the classically perfect actionlidO) . 

Finally we also tested our prediction regarding the degree of locality, and 
we found that indeed the exponential decay is much faster for the overlap 



c Also Dyy is ultralocal, but not an approximate GW fermion, so its change to Djv e i s rather 
drastic and the argument for a high level of locality does not apply to that case. 
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Figure 5: The "meson" dispersion relations for different types of overlap ferraions in the 
2- flavor Schwinger model at = 6. 

SO-HF than for the Neuberger fermion, see Fig. o and Table (jl. H 
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Figure 6: The level of locality for the Neuberger fermion compared to the overlap SO-HF: the 
phtcay of the free couplings in all directions (left) and the decay of the "maximal correlation" 
E3 over a distance r at (3 = 6 (right). The width of the "cones" in the left figure also shows 
how well rotation invariance is approximated. 



3 Overlap fermions and the doubling problem 

As we mentioned in Sec. 2, overlap fermions are free of doublers at least in a 
smooth gauge background. It should be clarified, however, that the overlap 
formula itself does not remove any doublers. They have to be removed before, 
in Do, by something like a Wilson term; in our ansatz Do — p^7^ + A, the 
scalar term A is crucial for that purpose. Then the overlap formula cures the 



d Rcf. J24 ] suggests to use non-standard operators (with (R) being subtracted from (D 1 }) 

even for the free fermion, which does improve e.g. the scaling of the free Neuberger fermion. 
However, one then deals with a non-local fermion, even though it is used only in an indirect 
way, which raises conceptual questions concerning the continuum limit. An exception ia-fche 
perfect fermion at R = 0, where for instance the axial anomaly is reproduced correctly Ej. 
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perfect 


overlap-HF 


Neuberger 


d = 2 


r P 


1.268 


1.255 


1.930 


d = 2 


r\ 


0.871 


0.888 


1.248 


d = A 


r P 


1.635 


1.519 


2.530 


d = 4 


r\ 


1.187 


1.109 


1.708 



Table 1: The characteristic radius r\ 
for various types of free GW fermions. 



(X] i J'M a 0l a;2 )/(X) : J'M a 0l) and (analogously) r p 



chiral symmetry again in the sense of the GWR, while doublers are supposed 
not to be re- introduced. 

We keep on referring to a point-like GW kernel R x ,y = ^-3x,y, where the 
mass parameter [i > can be chosen. It is the center of the circle through zero 
that the spectrum of the Dirac operator is mapped on by the overlap formula 
(which also involves the parameter f£) . Let us focus on the eigenvalues of Dq 
close to the real axis: the small (almost) real eigenvalues have to be mapped 
on (the vicinity of) zero, whereas the eigenvalues with large real parts have to 
be mapped on the opposite arc, i.e. on (the vicinity of) 2/i. For weak up to 
moderate coupling, there is a (statistically) safe interval where \i can be chosen 
such that this separation is achieved. Inside this interval, fi may be optimized 
with respect to various criteria; the optimal fi tends to rise with increasing 
coupling strength — since it has to adapLto the mass renormalization — as 
has been observed by optimizing localityEj or minimizing the mapping effect 
E3. For instance, in the Schwinger model such a safe interval still exists at 
j3 = 2, as Fig. g illustrates. 

But how about really strong coupling ? At some point, the eigenvalues of 
typical configurations spread all over inside a certain area for simple (short- 
ranged) operators Do. This is the case for the Wilson fcrmkin and for the HF 
in the Schwinger model at /3 = 1, and for QCD at (3 = 5,E3 as Fig. shows. 
Wherever we choose [i, mappings to the "wrong" arc will occur frequently; the 
possibility of "wrong" projections implies that both, the doubling problem (too 
many projections to the left arc) as well as mass renormalization (too many 
projections to the right arc) are back. Hence at really strong coupling, both 
of these problems return for the overlap fermions based on some simple Dq, 
as was also revealed by a strong coupling expansion cB H. Probably the only 
way out of this problem would be an excellent approximation to a (classically) 



e The question if also locality breaks down at some point of relatively slpang coupling, where 
the overlap formula is still applicable, is currently under investigation E3. For instance, one 
can check this for the Neuberger fermion in the Schwinger model at f3 < 2, or in QCD at 
Si 5.6. 
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Figure 7: Typical spectra of minimally gauged HFs and Wilson fermions (with the hopping 
parameter of the free fermion) in the Schwinger model at j3 = 1 (left) and in QCD at (3 = 5 
(right). This illustrates situations where the overlap projection is not applicable any more. 

perfect fermion, which is, however, very difficult to construct and implement. 

It is now of interest to get an idea where this transition occurs. Fig. g shows 
the Wilson spectra for typical (quenched) QCD configurations at /3 = 5.4 and 
(3 = 5.6. We see that the existence of a safe interval for /j, seems to set in in 
between. For the HF the situation is similar, just a little better. Hence 
j3 = 6 can be regarded a really safe regime with respect to this issue — Fig. || 
shows that an isolated left arc is present — and also with respect to locality 
ca, as we mentioned before. 

We therefore concentrate on QCD at (3 = 6 in the next Section, and we 
show that we can construct HFs, which approximate the GWR well at that 
coupling strength. This is the crucial ingredient needed to carry on the program 
of improved overlap fermions — described in Sec. 3 — to QCD. 

Note that in d = 4 the notorious square root has to be evaluated by some 
iterative method, and starting from a good GW fermion approximation instead 
of the Wilson fermion is also highly profitable with that respect: starting in 
the right vicinity speeds up the convergence a lot, as we .see for instance from 
the perturbative expansion of the square root around fi. 113 



' From our experience, even lattices as small as 4 4 provide a reliable insight into such ques- 
tions as the applicability of the overlap formula and the level of approximation to the GW 
circle. A disadvantage is that the left arc close to zero is missing — not due to the lattice 
action but just due to small volume. Therefore we also evaluated this subset of the spectra 
(the eigenvalues with smallest real parts) on 8 4 lattices, see below. 
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Figure 8: Typical spectra of the Wilson fermion in QCD on a 4 4 lattice at j3 = 5.4 (left) 
and j3 = 5.6 (right). We recognize the transition to the regime where the overlap projection 
is statistically safe (for 1.5</i<2j. 



4 Approximate Ginsparg- Wilson Fermions for QCD 

Our approach to construct a short-ranged approximate GW fermion for QCD is 
to stay with the couplings of the truncated perfect free fermion and gauge it by 
hand, using just very few new parameters to go beyond the "minimal gauging" . 
This concept was successful in d — 2, and in d = 4 we already know that the 
free HF is doing well in scaling, approximating the GWR and approximating 
rotational invariance (the latter is also checked at strong coupling), see Sec. 
1.3. Alternatively, one may tr y t o minimize the GWR violation directly within 
a limited set of parameters EZIE2I, or undertake a new effort to parameterize an 
(approximate) classically perfect actionEx 

We are confident that our free HF couplings already provide a good scaling, 
so the issue is to find a suitable gauging in the sense that the GWR violation 
is small. As our criterion, we compute the spectra on small lattices and try 
to arrange for them to be close to a GW unit circle (for typical quenched 
configurations at f3 — 6). 

As we see from Fig. |9|, the minimally gauged HF suffers from mass renor- 
malization almost as much as the Wilson fermion. On the other hand, the 
right arc is excellent (see Fig. [7]), but less important. Our first step beyond 
minimal gauging is the use of fat links: each link in a given configurations is 
substituted as 



link 



(l-a)Unk+- [ ^ 



staples 



{a e IR) 
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Figure 9: The left arc of typical spectra in QCD at (3 = 6 on an 8 4 lattice for the Wilson 
fermion and the minimally gauged HF. We show the 600 resp. 300 eigenvalues with smallest 
real parts, and we see that here the left arc is manifestly isolated. 

This is computationally cheap, and since we perform just one such fattening 
step we do not need to project back onto the gauge group (in contrast to Ref. 
[[32|])- The results for different a are shown in Fig. [To] (left). As we observed 
already in d — 2, a strongly negatpe a is required if one wants to remove the 
mass renormalization in this way 113. However, then the spectrum moves far 
away from the GW circle, so we do not recommend this way to approach the 
chiral limit. Positive a increases the pion mass further, but it makes the shape 
of the spectrum more circle-like. We are going to take advantage of that. 

Next we attach an amplification factor 1/u (u < 1) to each link to com- 
pensate the mean suppression by the gauge field. This is related in spirit to 
tadpole improvement Eil, but it can also be viewed as directly generalizing the 
tuning of the Wilson hopping parameter. For our value of (3 = 6, we reach 
criticality for the minimally gauged HF at u ~ 0.8, see Fig. [To] (right). Once u 
is fixed, its inclusion is computationally for free, and it does lead already to a 
decent approximation of the GW circle. 

In a next step, we include fat links with positive staples — a = 0.3 is a 
good value — and use again the critical link amplification parameter, which 
now amounts to u ~ 0.76. Indeed, this helps to move the eigenvalues closer to 
the GW circle, as Fig. [n] (left) shows. 

Still one would like to further reduce the imaginary part of the eigenvalues; 
in particular the upper (and lower) arc still calls for improvement. Considering 
the structure of Dhf(x, y, U) — p^(x, y, C/)7 M + A(x, y, U), we recognize that 
Pfj, is responsible for the imaginary part, so we multiply a damping factor v < 1 
on each link only in the vector term p^. The scalar term A, which controls 
the left and right arc already successfully, remains untouched. The optimal 
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Figure 10: The effect of fat links (left, u = 1) and link amplification (right, a = 0,u = 0.8j 
for a HF spectrum in QCD at /3 = 6 on a 4 4 lattice. 

value for the new parameter is v ~ 0.9 without fat links, and v ~ 0.92 at 
a = 0.3, always at critical u (which remains practically unchanged). Now 
also the upper arc follows the GW circle closely, and the fat link helps the 
eigenvalues to spread less around the circle, see Fig. [n] (right). This is the 
best approximation achieved so far, and we are confident that this is about the 
optimum that can be achieved with O(10) independent parameters. In fact, 
even if we include fat links we are using just 10 independent parameters. 

Of course one might still try further parameters, but they should only 
be included if they really lead to a significant progress. One could include 
terms with a new Dirac structure, and we are currently testing the clover 
term: its inclusion (with a positive coefficient) helps to improve the physically 
important arc near zero a little, but it distorts the opposite arc (in the sense 
that the eigenvalues fluctuate stronger around the right half-circle) . Generally 
the clover term tends to attract the eigenvalues closer to the real axis, which is 
also known from the Wilson fermiona. Hence the optimal value of v increases 
a little: for instance at clover coefficient 0.15 (and a — 0.3) it amounts to 
v ~ 0.94, while u ~ 0.767 is critical. 

The consideration of thp "magnetic mass" tub suggests that also a term 
oc 7^7iy7p could be useful t-3 (we refer to the Pauli term a B/2jjib in the low 
energy expansion). 

What is computationally simple and perhaps promising is an extension of 
the fat link to include also (selected) paths of length 5. All this is currently 
under investigation and the results will be reported in Ref. [p9|]- 
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Figure 11: QCD spectra on 4 lattices for the HF at critical link amplification parameter u: 
with and without fat links (left) and typical results if we further include the kinetic damping 
parameter v at its optimal value (right). We also show the "continuation" near in each 
case on an 8 4 lattice with the same parameters, everything at [3 = 6. 

5 Summary and Outlook 

Our concept is to first construct a short-ranged approximate GW fermion with 
a good scaling behavior and quasi rotational invariance. Inserting this fermion 
into the overlap formula makes the chirality exact. Since the alteration is small, 
the good scaling and approximate rotational symmetry is essentially preserved, 
so that we obtain an improved overlap fermion. For the same reason, such an 
overlap fermion has a high level of locality and its iterative evaluation is fast. 
This program was tested extensively in the Schwinger modelliJ. Furthermore, 
for very good approximate GW operators the interval (for the mass parameter 
/x) where the overlap formula can be used safely is extended, and such an 
interval exists up to stronger coupling compared to the Neuberger fermion 
(which inserts the Wilson fermion). ._. 

We add that the same mechanism applies to domain wall fexpjons E3: 
inserting an improved 4d fermion may induce the same advantagescO. 

For QCD, we use a truncated perfect free HF as our building block, and 
construct a good gauging in the sense of a small GWR violation. At j3 = 6 
this is achieved, as we observed from the spectrum: it is close to a GW circle, 
hence the fermion approximately obeys the GWR with a point-like kernel R. 

In our construction, we start form the minimally gauged HF, remove the 
additive mass renormalization by a critical amplification factor for each link, 
and suppress the vector term only in order to damp the imaginary part of 
the eigenvalues, so that the spectrum follows the shape of a GW unit circle. 
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Fat links help to reduce the fluctuations of the eigenvalues around that circle, 
and hence provide an approximation on a satisfactory level. It seems quite 
optimal for a set of about 10 parameters. Perhaps extended staple terms lead 
to some further progress, and the clover term helps a little with respect to 
the fine resolution close to zero, which we now focus on. In any case, a good 
approximation to a GW fermion in QCD has heen accomplished, and it is 
ready to be inserted into the overlap formula (|5|) lj. 
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